Initial spin periods of neutron stars in supernova remnants 
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ABSTRACT 

We present estimates of initial spin periods, Pq, for radio pulsars associated 
with supernova remnants. By using the published data on 30 objects, we were 
able to derive a reliable estimate for the initial spin period, assuming standard 
magneto-dipole spin-down (braking index n = 3), in many cases. Our set of 
estimates is still not sufficient to infer the exact shape of the initial period distri- 
bution. However, we show that a gaussian distribution with mean and deviation 
~ 0.1 s is consistent with our results, while flat, wide distributions and very 
narrow ones are disfavored. 

Subject headings: neutron stars 
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Introduction 



Studies of neutron star (NS) properties and evolution is a fruitful field of research 



uniting astrophysics and fundamental physics (see a genera 



example, in the introductor y part of the book 



Lattimer and Prakash 



Haensel et al. 



review and references, for 



2007 and a shorter review in 



2004 1 . In the map of this research area quantum chromodynamics 
is linked to observations of cooling NSs; general relativity with pulse profiles of X-ray 
sources; strong field electrodynamics with astrophysics of radiopulsars and magnetars, and 
so on. Despite significant progress, some parts of this picture are still obscure. One of 
them is related to NSs initial parameters. The knowledge of initial properties of compact 
objects is of crucial importance not only for understanding their evolution and modelling 
their observational appearance, but also because the parameters of newborn NSs carry the 
imprint of the still poorly understood supernova (SN) explosion mechanism. 

For a NS the main initial parameters are the mass (M), dipolar external magnetic field 
(B), spin period (P) and spatial velocity (v). In addition, the angle between the spin and 
magnetic axis, the initial configuration of the crustal magnetic field (including a toroidal 
component), and other parameters are important for the thermal and/or magneto-rotational 
evolution and for fixing the different observational manifestations of NSs. 

There is strong, ongoing effort by many research groups in the attempt to advance our 
knowledge of NSs initial parameter distributions. This can be done by different techniques. 
Direct observations, unfortunately, are impossible, but observational data can be used 
together with assumptions about the NS evolution to estimate a posteriori the initial 
parameters for individual ob jects. Another approach is related to population synthesis 



(IPopov and Prokhorov 



20071 ). In this case an evolutionary scenario is applied to follow 
numerically a large sample of sources, selection effects are taken into account, and results 
are confronted with observational distributions. Finally, initial parameters can be calculated 
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in the framework of some theoretical model (for example, models of SN explosion) 



Popul ation synthesis modelling wa s used to study the initial properties of radiopulsars 



(PSRs) by iFaucher-Giguere and Kaspil (120061 ). who considered several possible initial 



distributions o 



m 



Popov et al. 



NS pa rameters in the case of standard PSRs. Another approach was used 



( 120101 ). where in the framework of a decaying magnetic field, several types 
of isolated NSs were studied, with a focus on the properties of the initial magnetic field. A 
direct numer ical model is used , for example, to derive initial mass distribution of compact 



objects (e.g. 



Fryer et al. 



20111 ). Another exa mple is related to t he calculation of the kick 



(recoil) velocity of newborn compact objects (jScheck et al 



20061 ) 



In this note we address the problem of determining the initial spin period Pq of NSs. 
This is done by using data on the "true age" of NSs, as obtained by observations of 
associated SN remnants (SNRs) and applying a standard model of spin evolution. In the 
next section we present our sample and results. Then we provide a discussion of several 
issues related to our study, and, finally, conclude. 



2. Sample and initial periods 



We collected a sample of 30 NSs with well-measured P and P, associated with SNRs 
or/and pulsar wind nebulae (see Table 1 for a complete list). For all of them there are 
estimates of SNR ages based on different assumptio ns (in the case of the Crab n ebula and, 



probably, Gil. 2-0. 3 historical ages can be used, see 



Green and Stephenson 



and period derivatives are taken from the ATNF catalogue ( Manchester et al. 



2003). 



2005 



'eriod 



In 6 cases the independently estimated SNR age, tsnr, is consistent with the spin-down 



x The on-line catalogue is available at http://www.atnf.csiro.au/people/pulsar/psrcat/ 



-5 - 



age, r sc j = P/2P, calculated for zero initial spin periods and braking index n = 3 (n = vi)/v 2 , 
where v — 1/P is the spin frequency), and so no estimate of Pq is possible in the standard 
picture, except that it has to be Pq <C P. Note, that this does not formally mean that 
Po — 0. The initial period for these sources can range from a few ms up to tens of ms (as 
in, for example, PSR J1119-6127 and PSR J1846-0258). In plotting the estimated initial 
spin period for these sources (Fig. [Q we assume P < 0.1P. 

In three cases only upper or lower limits on a SNR age are available. For 15 SNRs age 
ranges are reported in the literature. For most of them it is possible to obtain a range of 
initial spin periods and for some the resulting ranges for P are narrow. In several cases r S( j 
is in between the upper and lower limit of tsnr- Correspondingly, only an upper limit on 
Po can be derived. Finally, in 6 cases the SNR ages are relatively well determined, and in 
one case the uncertainty is tiny in comparison with the spin-down age of the pulsar. This 
allows us to obtain firm estimates of Po for a given braking index. 

Initial spin periods are estimated from the usual expression 



Pn 



P' 



(n - 1)PP"- 2 tsnr 



l/(n-l) 



where P and P are the present values of the spin period and period derivative, and n is 
the braking index. For n = 3 the equation reduces to the usual magneto-dipole formula 



with B cx These are the values presented in Table 2 and Fig. 1. The field is 

assumed to be constant during the lifetime of a NS, as well as the angle between spin 
and magnetic axes. This assumption is well justified, as all sources are young and do not 
possess strong magnetic fields. So, it is very unlikely that field decay could play a major 
role in their evo l ution . Our estimates are in good correspondence with those provided by 



Migliazzo et al. 



(120021 ) . However, our list is significantly broader, although it is based on 



less detailed observational data. 
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2.1. Notes on individual objects 



N157B. The source is s i tuated in the Large Magellanic Cloud. The pulsar was 
discovered by 



Marshall et al 



1998) and the age of the associated Crab-like SNR is 



estimated in IWang and Gotthelj (119981 ) as ~5000 yrs basing on X-ray data (ROSAT and 
ASCA). Uncertainties in the age are not reported. The SNR age is in good correspondence 
with the characteristic age of the pulsar. 



G292.2-0.5. Both for free expansion and Sedov stage 



Pivovaroff et al 



(1200 ll ) conclude 



that the age of the nebula can be in correspondence with the radio pulsar characteristic age. 



GO. 9 +0.1. This is a comp osite SNR. Age estimates reported by lAharonian and et al. 



(120051 ) and 



Porquet et al. 



( 120031 ) are not precise, a few thousand years. This is compatible 



with the characteristic age of the pulsar. 

G359.23-0.82. This is a pulsar wind nebula. There are no reliable age estimates for 
this object which can be used for deriving the radio pulsar initial spin period. Formally, the 
age of the nebula is consistent with the pulsar characteristic age. 

Kes 75. Basing on the detailed analysis of the distance towards this object 



Leahy and Tian 
radio pulsar age 



( 120081 ) conclude that the age of the SNR is in correspondence with the 



G54. 1+0.3. This is a Crab-like nebula. ICamilo et al.l ( j2002a( ) reported that no 
independent reliable age estimate exists for this source. Roughly, the age of the nebula is 
consistent with the pulsar spin-down age. 



CTA 1. 



to 



Slane et al. 



Hafpern et al. 



( 120041 ) report an age of 13,000 yrs for the SNR age, referring 



(120041 ) . T his value correspon ds to a distance to the remnant of 1.4 kpc. The 



uncertainty, according to 



Slane et al. 



( 120041 ). is ±0.3 kpc, which translates in to ±2800 yrs 



for the uncertainty in the age. 
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3C 58. In a recent pape r 



to earlier studies. 



Slane et al.l (120081 ) report the value of 5400 yrs referring 



Bietenholzl (120061 ) gives 7000 yrs with a 3cr lower limit of 430 years 



Rud ie and Fesenl ( 120071 ) suggest a value of "about few thousand years" . Finally, 



Chevalier 



(120051 ) has 2400 ± 500 yrs. In our study we use the interval 4300-7000 yrs. 



S147. Old age estimates (see, fo r example 



Ng et al 



Kundu et al. 



1980)) suggested a SNR age 



( 120071 ) gave arguments in favour of shorter 



of ~ 80 - 200 kyrs. In a recent paper 
age, ~ 40 kyrs. The kinematic age of the pulsar is also close to this value. 



0540-693. Several estimates are available for this source. Manchester et al. (1993) give 



760 ±1 yrs referring to older studies, in particular to 



(ll989|). 



Reynolds ( 



762 ± 100 yrs. 



Reynolds! (119851) and 



Kirs 



1985 1) presents an esti mate of 800- 1 100 y ears. iKirshner et al. 



Williams et ai- 



mer et al. 



(120081 1 and 



Park et al. 



19891 ) give 



( 120101 ) provide wide lists of references 



and age estimates for this source. For our study we use an interval 660-1100 yrs. 



estimate for this object. They refer to 



Monogem Ring. iThorsett et al. l ( 2003) give a s i ngle y alue, 86,000 yrs, as the age 



Plucinsky et al. 



(119961 ) . where this value corresponds 



to a distance of 300 pc. However, possible distances range from 100 to 1300 pc. This can 
translates in to a 29,000-371,000 yrs interval for the age estimate. These authors state that 
small and large distances are not very probable, but d = 600 pc is certainly as good as 300 
pc. This distance estimate (600 pc) corresponds to the age 170,000 yrs. So, we use the 
interval 86-170 kyrs. 



Puppis A. iGotthelf and Halpernl ( 120091 ) refer to IWinkler et al.l (11988), where an 
est imate of 3700 ± 300 yr s is giv en. We use in our analysis the range 3400-4100 yrs reported 



by 



Gotthelf and Halpernl (120091 ) . 



Vela. For the Vela pulsar the upper limit for the age is taken from 



Aschenbach et al. 



( 119951 ). These authors used d = 500 pc, a value which is currently revised towards lower 
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values which result in a smaller upper limit for the age. The lower limit for this SNR is 
not very certain, and we use the value 11 kyr. Slightly smaller values are possible, but this 
would not change our conclusions about the initial period significantly. Strictly speaking, 
the age of the Vela remnant can be considered just in correspondence with the characteristic 
age of the pulsar. Still, we provide a full range for age estimates reported in the literature, 
and use the source in our analysis. 



G292.0+1.8. 



Gonzalez and Safi-Harbl (120031 ) give the range 2400-2850 yrs, which we 



use for our estimates. 



G320.4-1.2. The estimates of the SNR age given by lYatsu et al.l (120051 ) are based on 
standard assumptions about the SN explosion energy and ISM density. So, the age estimate 
is rather uncertain, probably the nebula's age is in correspondence with the pulsar age. We 
do not provide any estimated of Pq for this object. 



Roberts and Broganl ( 120081 ) is based 



G7.5-1.7. The estimate of 10 4 -10 5 yrs given by 
on numerical estimates for SNRs of this type, and as such is quite rough. The authors also 
suggest 50,000 yrs as a representative value. We use the full range above. 



G12. 8-0.0. We use the range of ages from 



Brogan et al 



(120051). In 



Dean and Hill 



( 120081 ) the properties of the radio pulsar are discussed in details, and arguments in favour 



of significant initial period (i.e., compatible with the present day one) are provided. 



G21.5-0.9 ISafi-Harb et al.l (120011 ) provide a list of references for age estimates of the 
SNR. We use the full reported range which is rather wide. 



W44. 



Harrus et al. 



(119971 ) present two possibilities. The first one gives a SNR age of 



5600- 7500 yrs with a central value 6500 yrs, the second one 19000-25000 yrs. In a recent 



paper 



Abdo and et al 



(120101 ) used an estimate ~ 20000 yrs. For our estimates we use the 



range 6.5-20 kyrs. 
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G65. 1+0.6. The large uncertainty in the age of the SNR is not influential on our 
results, as anyway the values are much smaller than the value of the characteristic age of 
the pulsar. 

CTB80. There is only an age estimate for the pulsar wind nebula. We consider this as 
a lower limit for the actual age of the pulsar, so only an upper limit for P Q can be obtained. 



G308. 8-0.1. ICaswell et al.l (Il992h report just an upper limit for r SNR because the 
value of the ISM density is uncertain (in particular it can be lower). This estimate is larger 
than r s d, and no conclusion about Pq can be reached. 



G106.6+2.9. 



Kothes et al. 



(120061 ) obtained an estimate for the pulsar wind nebula 



age. We use this as a lower limit for the age of the pulsar. 

Crab. The age of this SNR is well-known as it is related to a historical SN. 



G296. 5+10.0. IVasisht et al.l (119971 ) give an estimate of ~ 10 4 yrs for an exp l osion 



Kellett et al 



( 119871 ) 



energy E = 2 10 51 erg and an ISM density of 0.2 cm~ 3 . They refer to 
who provide an estimate 20,000 yrs for E = 1.5 10 51 erg and ISM density of 0.26 cm -3 . 
Both values are extremely small in comparison with the spin-down age of the pulsar. 



G315. 9-0.0. 



Camilo et al. 



(120091 ) calculate an age estimate basing on the assumption 
that the remnant is at the Sedov stage of expansion. Correspondingly, uncertainties are 
related to uncertainties in distance, ISM density, and SN explosion energy, which are 
however not discussed. 

G il. 2-0. 3. The ag e of this SNR is considered to be well established as it is a historical 



event ( jTorii et al.l 



19991 ). In addition, Sedo v phase estimates are consistent with the 



historical interpretation ( IVasisht et al. 



19961 ). 



Kes 79. Uncertainties in the age for this SNR are related to uncertainties in its 
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distance, which is of the order of ±1 kpc. Taking into account that the former are at 
the ~ 15 — 20% level and that the SNR age is orders of magnitude smaller than the 
characteristic age of the pulsar, we neglect the uncertainty. 



G78. 2+2.1. Exact uncertainty in the age is not reported by 
but it is anyway small and does not influence our results. 



Uchivama et all (120021 ). 



G114.3+0.3. Exact uncertainty in the age is not reported by lYar-Uyaniker et al. 



( 12004 ). but it is anyway small and does not influence our results. 



Table 1. Sample of PSRs associated with SNRs 



PSR 


SNR 


7"«jvr/10 3 yrs 


T sd /10 3 yrs 


Ref. 


J0537-6910 


N157B 


as the PSR 


4.9 


Wane and Gotthelf (1998) 


J1119-6127 


G292.2-0.5 


as the PSR 


1.6 


Pivovaroff et al. (2001) 


J1747-2809 


GO.9+0.1 


as the PSR 


5.3 


Aharonian and et al. (2005) 










Porauet et al. (2003) 


J1747-2958 


G359.23-0.82 


as the PSR 


25.5 


Camilo et al. (2002b) 


J1846-0258 


Kes75 


as the PSR 


0.73 


Leahv and Tian (2008) 


J1930+1852 


G54.1+0.3 


as the PSR 


2.9 


Camilo et al. (2002a) 


J0007+7303 


CTA 1 


10.2-15.8 


13.9 


Slane et al. (2004) 


J0205+6449 


3C58 


4.3-7 


5.4 


Slane et al. (2008) 


J0538+2817 


S147 


40-200 


618.1 


Anderson et al. (1996) 










Ne et al. (2007) 


B0540-69 


0540-693 


0.66-1.1 


1.67 


Williams et al. (2008) 


B0656+14 


Monogem Ring 


86-170 


110.9 


Thorsett et al. (2003) 


J0821-4300 


Puppis A 


3.3-4.1 


1489. 


Gotthelf and Habern (2009) 


B0833-45 


Vela 


11-27 


11.3 


Aschenbach et al. (1995) 


J1124-5916 


G292.0+1.8 


2.4-2.85 


2.85 


Gonzalez and Safi-Harb (2003" 


B1509-58 


G320.4-1.2 


6-20 


1.6 


Yatsu et al. (2005) 


J1809-2332 


G7.5-1.7 


10-100 


67.6 


Roberts and Broean (2008) 


J1813-1749 


G12.8-0.0 


0.285-2.5 


4.7 


Broean et al. (2005) 


J1833-1034 


G21.5-0.9 


0.8-40. 


4.9 


Safi-Harb et al. (2001) 
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Table 1 — Continued 



PSR 


SNR 


tsnr/10 3 yrs 


T sd /10 3 yrs 


Ref. 


B1853+01 


W44 


6.5-20 


OA O 

20.3 


Harms et al. (1997) 


J1957+2831 


Goo. 1+0. 6 


a r\ i /in 

40-140 


1568. 


lian and Leahy (2006) 




U 1 BoO 


> lo 


10 ( . 


Gastelletti et al. (zOOoj 


B1338-62 


G308 8-0 1 


< 32.5 


12.1 


Caswell et al (1992) 


J2229+6114 


G106.6+2.9 


> 3.9 


10.5 


Kothes et al. (2006) 


B0531+21 


Crab 


0.957 


1.24 


SteDhenson and Green (2002) 


J1210-5226 


G296. 5+10.0 


10-20 


101817. 


Vasisht et al. (1997) 


J1437-5959 


G315.9-0.0 


22 


114. 


Camilo et al. (2009) 


J1811-1925 


Gil. 2-0.3 


1.6 


23.2 


Torii et al. (1999) 


J1852+0040 


Kes79 


6 


191502. 


Sun et al. (2004) 


J2021+4026 


G78.2+2.1 


6.6 


76.9 


Uchivama et al. (2002) 


B2334+61 


G114.3+0.3 


7.7 


40.6 


Yar-Uvaniker et al. (2004) 



Table 2. Spin parameters of PSRs in the sample 



PSR P s P B/10 12 G P s P /P 



J0537-6910 0.016 

J1119-6127 0.408 

J1747-2809 0.052 

J1747-2958 0.099 

J1846-0258 0.326 

J1930+1852 0.137 

J0007+7303 0.316 
J0205+6449 0.066 
J0538+2817 0.143 
0.143 

B0540-69 0.05 
0.05 

B0656+14 0.385 
J0821-4300 0.113 
0.113 

B0833-45 0.089 
J1124-5916 0.135 
0.135 

J1210-5226 0.424 
B1509-58 0.151 



5.18E-14 0.92 

4.02E-12 41. 

1.56E-13 2.9 

6.13E-14 2.5 

7.08E-12 48.6 

7.51E-13 10.3 

3.6E-13 10.8 

1.94E-13 3.6 

3.67E-15 0.73 

3.67E-15 0.73 

4.79E-13 5.0 

4.79E-13 5.0 

5.5E-14 4.7 

1.2E-15 0.37 

1.2E-15 0.37 

1.25E-13 3.4 

7.53E-13 10.2 

7.53E-13 10.2 

6.6E-17 0.17 

1.54E-12 15.4 





~ 


<.p 


~ 


<.p 


~ 


< p 


~ 


<.p 


~ 


<.p 


~ 


< 0.163 


< 0.52 


< 0.029 


< 0.45 


< 0.134 


< 0.93 


> 0.118 


> 0.82 


< 0.039 


< 0.78 


> 0.03 


> 0.59 


< 0.183 


< 0.48 


< 0.113 


- 1 


> 0.113 


- 1 


< 0.016 


< 0.2 


< 0.054 


< 0.40 


> 0.004 


> 0.03 


0.424 


- 1 
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Table 2 — Continued 



PSR P s P P/10 12 G P s P /P 



J1809-2332 


0.147 


3 


,44E- 


•14 


2.3 


< 0.136 


< 0.92 


J1813-1749 


0.045 


1 


.5E- 


13 


2.6 


< 0.043 


< 0.97 




0.045 


1 


.5E- 


13 


2.6 


> 0.031 


> 0.69 


J1833-1034 


0.062 


2 


.02E- 


13 


3.6 


< 0.057 


< 0.91 


B1853+01 


0.267 


2 


.08E- 


13 


7.5 


< 0.221 


< 0.83 




0.267 


2. 


.08E- 


13 


7.5 


> 0.036 


> 0.14 


J1957+2831 


0.308 


3. 


.1 IE- 


■15 


0.99 


< 0.3 


< 0.99 




0.308 


3 


HE- 


■15 


0.99 


> 0.29 


> 0.95 


B1951+32 


0.04 


5 


.84E- 


15 


0.49 


< 0.036 


< 0.91 


B1338-62 


0.193 


2 


.53E- 


13 


7.1 






J2229+6114 


0.052 


7. 


.83E- 


■14 


2.0 


< 0.041 


< 0.79 


B0531+21 


0.033 


4. 


.23E- 


13 


3.8 


0.016 


0.48 


J1437-5959 


0.062 


8. 


.59E- 


15 


0.74 


0.055 


0.9 


J1811-1925 


0.065 


4. 


.40E- 


■14 


1.7 


0.062 


0.97 


J1852+0040 


0.105 


8. 


.68E- 


■18 


0.03 


0.105 


- 1 


J2021+4026 


0.265 


5. 


.47E- 


■14 


3.9 


0.254 


0.96 


B2334+61 


0.495 


1 


.93E- 


13 


9.9 


0.45 


0.91 
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3. Discussion 

For our estimates we use a braking index n — 3. Additional checks were made for 
other values, 2 < n < 10, and we found that the derived values of Pq are not significantly 
changed in many cases. For n = 2 P may be larger by a factor up to 2, but these changes 
do not modify the general picture shown in Fig. 1 significantly. This is an expected result, 
since the obje cts under consideration are very yo ung, and is in correspondence with the 



conclusions by iFaucher-Giguere and Kaspi 



(120061 ) that varying n around the standard value 



has a weak impact on the parameters of the pulsar population. However, with growing n 
more and more SNR ages are found to be in contradiction with the spin-down ages. Already 
for n = 4 in several cases r s d become smaller than tsnr- 



Recently, an e-print by IZhang and Xid (120111 ) appeared in which a correlation between 



magnetic field and the ratio of the "true" (i.e. tsnr) and spin-down ages was proposed. 
The authors explain this correlation as due to field decay. We present a similar plot, i.e. B 
vs. ages ratio, in Fig. 2 using data from Table 2. However, we believe tha t the existence o f 



such a correlation is an artifact, and the reason for its appearance is that 



Zhang and Xie 



(120 111 ) neglect the initial spin period. This conclusion can be illustrated by plotting B vs. 
the difference P — P . Clearly, in the sample of young objects (associated with SNRs) those 
with smaller magnetic fields have initial spin periods closer to the present day periods. 

Naturally, objects with lower field, like the "anti-magnetars" associated to some central 
compact objects in SNRs, have present spin periods very close to the initial values. Then, 
the usual estimate of the spin-down age assuming P ~ is invalid, as the true age is very 
different. 

We cannot exclude that there was some field evolution in the objects under discussion 
(note that field evolution is undistinguishable from the evolution of the angle between 
spin and magnetic axis), but we see no evidence for it, contrary to the conclusion by 
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£ C-&< □ □ 



□ 



<- 



-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

log Po 



Fig. 1. — Magnetic fields and initial spin periods for 22 PSRs (squares). In several cases only 
upper limits are available (shown with arrows). For the group of six pulsars with SNR ages 
in correspondence with spin-down ages we use different symbols (crosses) and assume that 
Po < 0.1P, where P is the present period. All values are estimated for standard magneto- 
dipole losses in the case of constant field and angle between spin and magnetic axis. 
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Fig. 2. — A plot similar to the one in 
"true" age (t snr ) and spin-down age. 
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Not surprisingly, the assumptions made by 



( jZhang and Xie 



Zhang and Xid (1201 if ) led them to suggest 



2012]) that field decay is so strong and important that in their youth most 
of standard pulsars have magnetar properties, B ~ 5 10 14 G at the age ~ 0.4 yrs. We 
believe that this is not needed if non-zero initial spin periods are accounted for. 



Oscillations of braking indices were also proposed by iBiryukov et al.l ( 120121 ) on a time 
scale ~ 10 3 - 10 4 yrs. If this is valid also for young pulsars, then our estimates of Pq must 
be corrected. But still, non-zero Pq will remain the main reason to explain discrepancy 
between spin-down ages and SNR ages. 

In Fig. 1 we see that under the assumptions we made no correlation between P Q and B 
is visible. This is in correspondence with standard assumptions made in different kinds of 
models. We did not try to figure out a possible shape of the initial spin period distribution 
(for example, if it is a multi-peak distribution or not — in the future, probably, it will 
appear that, for example, so-called "anti-magnetars" form a separate subpopulation, and 
so they do not fit the same single-mode distribution as normal radio pulsars; if gaussian is 
better than log-gaussian, symmetric distributions better than asymmetric, etc.). However, 
the data we obtained can be used to check assumptions about initial spin distributions. For 
example, obviously very narrow distributions are not in correspondence with the data as 
there are pulsars with Pq from milliseconds to hundreds of milliseconds. 

In order to derive information on the distribution of the initial periods from the values 
derived from observations, we run a number of Kolmogorov-Smirnov (KS) tests, using 
different forms of the distribution function. In particular, we tested gaussian distributions 
with different parameters and a top-hat (i.e. constant in a range and zero outside) 
distribution. For this calculations we used data on 11 PSRs. Six of the are those with an 
exactly determined Pq, to which we added five objects with small uncertainty for the initial 



-20 - 



period value using average SNR age. 

Two representative examples are shown in Figs. [3H33 The gaussian distribution of Fig. 
[3] gives a KS significance level of ~ 0.66, indicating that the data and model have quite 
consistent cumulative distributions. On the other hand, the top-hat distribution has a very 
low value of the KS significance, < 0.01, and should be rejected. 



4. Conclusions 

In this paper we discussed a sample of radio pulsars associated with supernova 
remnants. This association allows to obtain estimates of initial spin periods in > 20 cases. 
The obtained distribution is rather flat in wide range from ~ 10 msec to hundreds of msec. 

We thank Andrea Possenti for some helpful comments on the paper. SP also thanks 
Anton Biryukov and Andrei Igoshev for discussions. The work of S.P. was supported by 
RFBR grants (10-02-00599, 12-02-00186) and by the Federal programm for scientific and 
educational staff (02.740.11.0575). RT acknowledges financial support by INAF through a 
PRIN 2011 grant. SP thanks the University of Padova for hospitality. 
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Fig. 3. — The data cumulative distribution function (full histogram) compared with that of 
a gaussian distribution with average 0.1 s and standard deviation of 0.1 s (dotted line). The 
vertical line shows the value of P where the difference is largest. 
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Fig. 4. — Same as in Fig. |3]for a top-hat distribution in the range 0.001 s < Pq < 0.5 s. 
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